Extremely fast focal-plane wavefront sensing for extreme 

adaptive optics 

Christoph U. Keller a , Visa Korkiakoski a , Niek Doelman 6 , Rufus Fraanje c , Raluca Andrei c and 

Michel Verhaegen c 

a Leiden Observatory, Niels Bohrweg 2, 2333CA Leiden, The Netherlands; 
6 TNO Science and Industry, Stieltjesweg 1, 2628CK Delft, The Netherlands; 
c Delft Center for Systems and Control, Mekelweg 2, 2628CD Delft, The Netherlands 

ABSTRACT 

We present a promising approach to the extremely fast sensing and correction of small wavefront errors in 
adaptive optics systems. As our algorithm's computational complexity is roughly proportional to the number 
of actuators, it is particularly suitable to systems with 10,000 to 100,000 actuators. Our approach is based on 
sequential phase diversity and simple relations between the point-spread function and the wavefront error in the 
case of small aberrations. The particular choice of phase diversity, introduced by the deformable mirror itself, 
minimizes the wavefront error as well as the computational complexity. The method is well suited for high- 
contrast astronomical imaging of point sources such as the direct detection and characterization of exoplanets 
around stars, and it works even in the presence of a coronagraph that suppresses the diffraction pattern. The 
accompanying paper in these proceedings by Korkiakoski et al. describes the performance of the algorithm using 
numerical simulations and laboratory tests. 
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1. INTRODUCTION 

The computational requirements for extreme adaptive optics for the next generation of 30 to 40-meter class 
astronomical telescopes present a formidable problem. For an adaptive optics system with a phase corrector 
that has N degrees of freedom, classical approaches require a multiplication of a vector of length N with a 
matrix of size N x N, leading to a computational complexity of order iV 2 . For extreme adaptive optics, this 
vector-matrix multiplication is the major bottleneck, not only in terms of computational power but also in terms 
of memory requirements as both scale with TV 2 . Indeed, a system with N = 40,000 as envisaged for the EPICS 
instrument at the E-ELT, the classical reconstruction matrix would require about 6GB of memory (assuming 
4-byte single-precision floating point numbers). This data needs to be moved into the processing units at a rate 
of a few thousand times per second, making this issue alone a challenging data flow problem. 

Even when using parallel processing and distributed memory techniques, the most modern hardware imple- 
mentation would struggle to achieve the computational complexity and speed. Depending on the type of wave- 
front sensor and deformable mirror, sparse-matrix approaches,^ sequential one-dimensional reconstructions^ and 
Fourier reconstruction techniques^ can reduce the required processing power. Typically, these approaches work 
with conventional pupil-based wavefront sensors such as Shack-Hartmann and pyramid wavefront sensors. 

Here we present a focal-plane sensing algorithm with a complexity that is proportional to N log N in terms of 
required computing power and proportional to N in terms of memory requirements. It is based on 1) the close 
relationship between the PSF and the wavefront aberration in case that the latter is smalPH^, 2) sequential phase- 
diversitjdEl where the deformable mirror itself introduces the phase diversity and 3) a choice of the introduced 
phase diversity that minimizes the wavefront error and fortunately also minimizes the computational effort. 
Indeed, our algorithm only requires a few floating point calculations per wavefront resolution element and a 
single, two-dimensional, complex Fourier transform per update cycle. Our work combines knowledge from the 
areas of AO-corrected speckle studies^ and the phase-diversity efforts of GonsalvesGDH 
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We begin in Sect. [2] by calculating the monochromatic PSF of a weakly aberrated system to second order, 
which allows us to derive a correction to the first-order approximation that makes it significantly better. We 
continue in Sect.[3]by deriving the basic sequential phase diversity technique for our approximation. Sect.[4]shows 
how the computational complexity can be dramatically reduced from four Fourier transforms to a single one. 
The advantages of our approach are discussed in Sect. [5] the limitations are listed in Sect. [6j and applications 
and extensions are summarized in Sect. [7J In a separate paper in these proceedings, we have also studied in 
detail one version of our algorithm numerically as well as in the laboratory with a 37-actuator systemP 



2. MONOCHROMATIC PSF DUE TO WEAK ABERRATIONS 

The overall approach is based on the work by GonsalvesP who presented an analytical solution to retrieving the 
wavefront from two simultaneous images of a point source where one is in focus and the other one has a phase 
diversity, typically defocus, applied. Gonsalves^HU also introduced sequential phase diversity (SPD) where the 
phase-diverse images are acquired sequentially instead of simultaneously. We expand Gonsalves' work here to 
second order with a notation that is more closely related to the work that has been done for understanding the 
PSF of AO-corrected point sources.^ This expansion leads to a much deeper understanding of the approach 
suggested by GonsalvesP and implies a modification that significantly expands the range of the weak aberration 
approximation as shown in the accompanying paper by KorkiakoskiP 

2.1 First-order approximation of wavefront 

We start with the (scalar) electrical field of an electromagnetic wave in the pupil plane of an optical system, 

E(u, v) = A(u, v)e l * M (1) 

where A(u, v) is the amplitude and ^(u, v) is the phase in radians, both being functions of the spatial coordinates 
u and v in the pupil plane. The origin of the (it, v) coordinate system is the center of the aperture. Both A(u, v) 
and &(u,v) are real. Furthermore, as the absolute value of the phase (piston term) is not relevant, we can set 
the average of <& over the aperture to zero, i.e. 

J A(u, v)dudv = 0. (2) 

Note that A(u, v) does not have to be limited to values of or 1 but can have values in between such as used in 
apodized-pupil coronagraphsP-S 

As we do not rely on absolute intensity measurements and assume that we only deal with pure phase aber- 
rations, the amplitude of the electrical field is normalized such that 

J ' \E{u,v)\ 2 dudv = J \A{u,v)\ 2 dudv = 1. (3) 

Furthermore, we assume that the amplitude is even, i.e. 

A(u,v) = A(-u, -v) . (4) 

This is an assumption that will greatly simplify the subsequent calculations without being overly restrictive. An 
even amplitude or aperture function A(u,v) still describes most telescope apertures but also holds for shaped 
pupils and apodized pupils that are used to suppress diffraction effects around a point source.^ 

If the aberrations are very small, i.e. $ < 1, E(u,v) in Eq.Q can be approximated to first order by 

E(u, v) « E'(u, v) = A(u, v) [1 + v)] . (5) 



We note that the integral over the absolute value squared of E' is not unity anymore since 

\E'\ 2 dudv= / A(u, vf [l + v) 2 ] dudv = l+ / [A{u, v)$(u, v)} 2 du dv . (6) 



The error made in the normalization of the amplitude by this first-order approximation is therefore given by 
the variance of the aperture- weighted wavefront, A$ because of Eq.Q. This violation of energy conservation is 
the fundamental reason that algorithms based on Gonsalves^ first-order approximation are limited to extremely 
small phase aberrations and that the Strehl ratio of the corresponding PSF is always 1, independent of the 
wavefront variance, as will be shown below. 



2.2 Second-order approximation of wavefront 
Expanding Eq.([I]) to second order, we obtain 

E(u,v) w E"(u,v) = A(u,v) 

The integral 

\E"{u : v)\ 2 dudv= I A(u,v) 



1 + v) - ^(u, v) 2 



(7) 



dudv^l (8) 



l-^$(ti, v) 2 j +<&(u,vf 

up to second order for $ <C 1 since (1 — |) 2 w 1 — x) for x <C 1. Therefore, expanding Eq.jl| up to second order 
in the wavefront aberration $ will guarantee energy conservation up to second order. 

In the first-order approximation of the complex electrical field in the aperture, the real term is simply 1. 
In the second-order approximation, the real term of the series becomes 1 — |$(w,f) 2 , which will always be 
smaller than 1. The first-order expansion will therefore always overestimate the real term. We can correct for 
this overestimation in a statistical sense by replacing the general l$(it,w) 2 term in the complex field amplitude 
approximation with its average (over the pupil), 

11/" 1 

-$(u, v) 2 « - / A(u, v) 2 <S>(u, v) 2 du dv = -<j% , (9) 

where a\ is the wavefront variance in radians squared. We can therefore mitigate the problem in the first-order 
approximation of the pupil electric field by scaling the real part with 1 — <j%/2. The imaginary part will not 
require any scaling. 

We note that the square of this scaling factor, to second order in $ 

(i-4/2) 2 «i-4 (io) 

is the same as the Strehl ratio for small aberrations, the extended Marechal approximation.^ At the end of 
this subsection, we will provide an alternative explanation for this scaling factor that is closely linked to the 
second-order expansion of the wavefront aberration and provides a clear connection between this term and the 
Strehl ratio. 

2.3 Fourier transforms and symmetries 

Following Gonsalves^ we split the phase $ into odd and even terms such that 

$(u,v) = <S> (u,v) + $ e (u,v) (11) 

and 

$ {U,V) = -$„(-«,-«), $ e (u,«) =$ e (-U,-«) . (12) 

The corresponding approximate complex electric field amplitude in the focal plane is given by the Fourier 
transform of E"(u,v) in Eq.([7]). The coordinate system in the focal plane has its origin at the center of the 
perfect PSF without aberrations. We will not explicitly write the dependence on the focal-plane coordinates to 
make the equations more readable. Hence, the electrical field in the focal plane to second order in the phase 
aberration $ is 

e" = a + ia* (<j> e + <f> ) - -a * (<p e + tj> ) * (<p e + <f> ) , (13) 

where lower-case symbols are the Fourier transforms of the upper-case symbols, and * is the convolution operator. 

The Fourier transform of a real, even function is real and even, and the Fourier transform of a real, odd 
function is a purely imaginary, odd function. Therefore, <j) e and a are real and even, and <p a is imaginary and 
odd. 



2.4 Point-Spread Function 

The point-spread function (PSF) to second order in the aberration, p, is given by the absolute value squared of 
the complex field amplitude in the focal plane, 

p = e"-e"* (14) 
= a 2 (15) 
+ia [a * [4> e + <j) ) - a * (</>* + 0*)] 
+ [a * {<f> e + cj> )] [a * {4>* e + 4>* )] 

-~a[a* (<f> e + <f> ) * (<f> e + <f> ) + a * (ft + <f>* ) * + 0*)] , 

where * indicates the complex conjugate and we have dropped terms proportional to <£ 3 and higher. With A(u, v) 
being real and even, a = a* is real and a 2 is simply the PSF without aberrations. Since </>* = <p e and 0* = — (j) , 



Eq.(16) simplifies to 



p(x,y) = a 2 (16) 
+2ia(a * <j) ) 
+(a * <j) e ) 2 - (a * 4> Q ) 2 
-a (a * (f) * (f> + a * (f> e * cp e ) . 

This is the PSF to second order in the aberrations Note that the last term is mostly missing in publications 
that use the first-order approximation in the wavefront to derive an approximate PSF. 

Using the notation from Gonsalves^ work on phase-diversity in the weak aberration limit, 

v = a * 4> e (17) 
y = ia*(j) (18) 

and both v and y being real quantities, we can rewrite EqjlT] as 

p = a 2 + 2ay + y 2 + v 2 — a (a * <j) e * (p e + a * <f> * <j> ) . (19) 



Following again GonsalvesP we separate the PSF, p, into its odd and even parts, 

Po = 2ay (20) 
Pe = a 2 + y 2 + v 2 - a (a * 4> e * 4> e + a * 4> * (f> ) . (21) 

The last term for the even part of the wavefront is difficult to include in analytical solutions of the wavefront 
sensing problem as it includes convolutions of wavefronts with themselves and the aperture. However, it is the 
only term that can be negative and that is not zero at the center of the unaberrated PSF, the coordinate origin 
in the focal plane. The latter is due to our previous assumption that the phase, averaged over the aperture, is 
zero. It is indeed this term that leads to a reduction of the Strehl ratio as has already been pointed out by Perrin 
etalP 

Due to the requirement for energy conservation, the integral over the PSF has to remain constant, independent 
of the wavefront aberration. Therefore 

J y 2 + v 2 — a{a* (j) e * 4> e + a* <j) * <j> ) dx dy — . (22) 

Since the first two terms are always positive, the third term has to be largely negative. Neglecting this term 
therefore leads to a consistent bias in the estimate for the amplitude of v, which limits the applicability of 
Gonsalves' algorithm to very small phase aberrations. 



To simplify matters, we approximate Q(u, v) 2 with its average a%. <f>* <j> therefore becomes a delta function 
and the last term in Eq.(21 1 becomes 



a (a * <fi e * cj) e + a * <f) * <f> ) a\a 2 . (23) 
Finally, the equation for the even part of the PSF reduces to 

p e = (1 - a%)a 2 + y 2 + v 2 . (24) 

The even part is a modified version of the equation given by GonsalveiP where we correct the unaberrated 
PSF, a 2 , with the factor I — a 2 ,. Even though we do not know the variance of the wavefront aberrations in 
advance of the actual wavefront sensing, we can easily determine this correction factor as it is the same as 
the Strehl ratio for small aberrations in the extended Marechal approximation.^ We can therefore normalize 
the unaberrated PSF with the (observed) Strehl ratio, which is given by the maximum of p a bs- An equivalent 
estimate of the correction factor can be obtained from the requirement for energy conservation: a\ = J y 2 + v 2 . 
This correction to the equation given by Gonsalves^ significantly extends the amplitude of the aberrations over 
which the algorithm can be applied. Indeed, an rms aberration of about 1.5 radians can be toleratedP 

2.5 Odd part of wavefront error 

Once more following GonsalvesPy is easy to calculate from the odd part of the PSF, Eq.(21) via 



P° ~ 1 P° a mk\ 
la 2 a + e 

where e > avoids noise amplifications in places where a is very small. The odd part of the wavefront error, 
$o(u, v) multiplied with the aperture function A(u,v) is then obtained from the imaginary part of the inverse 
Fourier transform of y, 

A$ = 3Y~ . (26) 

2.6 Even part of wavefront error 

The magnitude of v, the Fourier transform of the even part of the phase function multiplied with the aperture 



function can be easily determined from Eq. ( 24 1 , 



\p e - (I - a 2 )a 2 - y 2 \ . (27) 

We use the absolute value here as this is only approximately correct, and therefore the argument of the square 
root has, in practice, also negative values. Preliminary experiments have shown this to be better than to set all 
negative elements to zero. The sign of v cannot be determined from a single PSF measurement. A convenient 
way to determine the sign of v is through phase diversitySEH an additional, known aberration is added, very 
often a defocus term, to figure out the sign of v. This is the topic of the next section. Once a sign has been 
determined, an inverse Fourier transform of v leads to the even part of the wavefront error multiplied with the 
aperture function, 

A$ e = $V . (28) 

Note that the inverse Fourier transform of v should be limited to the support of the aperture A. Any error 
in the choice of the sign may lead to Fourier components outside of the aperture, A (see Fig.[T]). Therefore, the 
sign should be chosen so that the inverse Fourier transform is limited to the aperture. 




Figure 1. The (signed) square root of the even part of the PSF due to the even wavefront error, v, (top row), along with its 
Fourier transform, V , (bottom row), which is an estimate of the even part of the wavefront error, "l>e. The three columns 
differ in terms of the accuracy of the sign of v: the correct sign (left column), all positive signs (middle column), and 20% 
wrong signs (right column). Note that sign erros lead to Fourier amplitudes outside of the aperture and a corresponding 
reduction of the amplitudes inside the aperture. 



3. SEQUENTIAL PHASE DIVERSITY 

3.1 Phase Diversity and Adaptive Optics 

Phase diversitjff^^ is a convenient way to determine the phase aberrations of images. Typically, two images 
are recorded simultaneously, often an in-focus and an out-of-focus image. Another approach, suggested by 
GonsalveSS consists in using a deformable mirror in an adaptive optics system to sequentially introduce an 
additional phase aberration. It is indeed easy to see how an adaptive optics system in normal operations can be 
used as a sequential phase diversity wavefront sensor. Imagine that at time tk we observe a partially corrected 
image that has an (unknown) phase aberration If we can assume that the deformable mirror can change 
its shape faster than the typical changes in the incoming wavefront that is being corrected, then the deformable 
mirror of the AO system in the previous update cycle had introduced an additional phase aberration of —A<f>k-i 
where A^^-i is the correction that was made by the deformable mirror in the last iteration. 

To use phase-diversity as a wavefront sensor in an adaptive optics system, one therefore has to know the 
previous and the current focal-plane images along with the phase change introduced by the deformable mirror 
during the previous update cycle. 

3.2 Choice of diversity 

In the weak-aberration regime developed above, phase diversity only has to help us find the sign of v. In most 
cases, one chooses a quadratic diversity as it is easily implemented by defocusing the imaging camera. When 



using a deformable mirror, much more elaborate choices of diversity can be chosen. Indeed, our goal is to find a 
diversity for every update cycle that minimizes the computational effort of determining the sign of v. 

If one uses a small diversity, it is easy to show^ that the even part of the PSF is 

p' e = (l-a'i)a 2 +y 2 + (v + z) 2 , (29) 

where z — a * 5<p e , where S(j) e is the even part of 5<j), which in turn is the Fourier transform of A$. If we know 
the introduced diversity, we can easily calculate z. Note that the introduced diversity does not have to be even, 
as is almost always assumed in phase-diversity work, since odd and even wavefront errors are nicely separated 
in the odd and even parts of the PSF. It is therefore sufficient to look at the even part of the introduced phase 
diversity. 

A simple approach to minimize the computational effort is the following: start by choosing an arbitrary sign for 
v and call it Vk-i- In the next update cycle of the AO system, we make the assumption that z — —v^-\ (note that 
this is not true since z is related to the introduced diversity, not what we think that we introduced. Therefore, if 
the choice of sign was correct, (vk — v/c-i) 2 = 0, and if the sign was wrong, (ufc — Vfc-i) 2 « Av\_ x . The prescription 
is then to monitor p e at every point in the focal plane. If the magnitude of (p e ,k ~~ Uk)( x ' v) > (Pe.fe-i — Vk-i)i x i u) 
then one flips the sign at location (x,y), i.e. sign Vk = —sign i>k—i- The problem with this simple approach is 
two-fold: 1) the algorithm in no way takes into account that the inverse Fourier transform of v should have a 
support limited by A; 2) once the correct sign choice has been made, the wavefront error will not be zero because 
of all the approximations one has made, and in particular the Fourier issue mentioned under 1). The choice in 
sign has again to be a random choice, which reduces the convergence speed by about a factor of 2 from what we 
would like to achieve. 

A significantly better approach is modeled after the one introduced by Gonsalves^ where we take into account 
that we also update the odd wavefront, and not just the even wavefront. The sign of Vk in iteration k is then 
determined by 

signt-fe = sign( W 2 _ 1 - v 2 - z 2 )/(2z) . (30) 

This approach has two significant advantages, possibly not realized by GonsalvesP 1) it deals with the Fourier 
support issue as Vk-i does not necessarily obey the Fourier support contraint, but z does, and 2) even if the 
correct sign choice was made for Vk-i, there will be a residual error, and the above equation also tries to estimate 
the best choice in this iteration based on taking into account the difference between Vk-i and z. 

The problem is equivalent to being able to finding the minimum of y — x 2 in the case that we can change x, 
but only measure y. We measure an initial value of x\ . Without any additional information, we must randomly 
choose the sign of the step that we are going to take to minimize y. Without loss of generality, we choose 
Sxi — — yx^, i.e. we assume that we are on the right side of x — 0. We then measure a value of x\. Obviously, 
if x\ > x\ , we chose the wrong sign. If we chose the correct sign, then the next step is 8x2 = iy^cf. But we 
would also like to know on which side of x = we are after having taken the step 6x% so that we do not have to 
again guess the sign in the next iteration. The correct sign for 6x2 can be estimated form the following relation: 

x \ = (%2 + 5xi) 2 = x\ + 2x25x\ + 5x\ . (31) 

Therefore the sign of X2 is given by the following equation 

x 2 = {x\ -x\- 5x 2 1 )/(25x 1 ) . (32) 

If X2 is positive, 5x2 must be negative, and vice versa. 

Even if the commanded value of 5x\ and the actually achieved step in x are not the same, the above equation 



is still correct when replacing 5x\ with the actual step z, which is then equivalent to Eq.(30|. This is particularly 
important in our case as we will always have errors in the choice of sign for some parts of v. Because of these 
errors, the deformable mirror will not achieve the correct shape, and the actual change will lead to a step size z 
that is different from what we commanded initially. 
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Figure 2. Example of the convergence of our algorithm in terms of the PSF (top), total wavefront error (second from 
top), odd wavefront error (second from bottom) and even wavefront error (bottom) along with the iteration number at 
the very bottom and the Strehl ratio at the top. 

4. REDUCING COMPUTATIONAL COMPLEXITY 

An obvious way to minimize the computing requirements consists in reducing the number of Fourier transforms. 
Due to the symmetry properties of Fourier transforms, the inverse Fourier transforms of y and v can easily be 
combined into a single, complex Fourier transform according to 

A(u, «)'$'(«, v) = ^(V + Y)+ 9f(V + Y) . (33) 

The calculation of z also requires two Fourier transforms to implement the pupil plane support restriction. 
Realizing that this 1) low-pass filtering in the Fourier domain is nothing more than a two-dimensional filter in 
the focal plane (a sine filter in the case of a top-hat aperture function A(u, v)) and 2) that this filter in the focal 
plane is strongly concentrated around the origin, we can perform this filtering with a brute-force convolution 
with a very limited kernel of about 5 by 5 pixels. This will add about 25 multiply-add calculations per actuator, 
which is negligible compared to the NlogN complexity of a fast Fourier transform. 

Since the limited support of V in the Fourier domain is nothing more than a low-pass filter applied to v, 
one can compare the estimate of v one has obtained and compare it to the convolution of v with the Fourier 
transform of A, which is v * a. If v and v * a have the same sign, then the chosen sign for v is correct, otherwise 
the opposite sign is a better choice. This post-processing step can again be implemented with a direct convolution 
with a very small kernel. 

Figure [2] shows 10 iterations of our algorithm with the direct convolution approximation for z and the post- 
processing of the sign choice as described above for a static, Kolmogorov-spectrum aberration with a starting 
Strehl ratio of 0.35. 

5. ADVANTAGES 

The obvious advantage of our approach is that it only requires a single, complex Fourier transform and a 
number of operations proportional to the number of actuators. In contrast to other Fourier reconstructorsP 
our approach needs no particular treatment of the edges as we are directly sensing the phase, not its spatial 



derivative. Computationally, it is therefore very fast and can easily cope with systems that have 10,000 to 100,000 
actuators. 

Some phase retrieval algorithms have a tendency to stall. This is not the case here as any increase in Strehl 
ratio will also mean that the weak-aberration assumption is more accurate, and therefore the algorithm will 
perform better. 

Our algorithm is also ideally suited for sensing a wavefront after an apodizing pupil coronagraph. The odd 
term and the second term in the even part of the PSF are both products with a, the complex field corresponding 
to the perfect PSF. As the main task of a coronagraph is to reduce the magnitude of a, these two terms will also 
decrease correspondingly. The estimation of the even part of the wavefront therefore becomes much easier since 
the second term contributes much less. 

We can also limit the radial extension of the PSFs that we work with, and can thereby limit the number of 
Fourier modes that are reconstructed. Even an obscuration from a coronagraph image mask can be taken into 
account if the corresponding low-order aberrations are well corrected. 

Finally, our approach easily implements any Fourier-based technique to make predictive corrections, e.g. 
shifting the entire wavefront in one direction, or take into account the influence function of actuators as long as 
it can be described by a convolution. 

6. LIMITATIONS 

Our wavefront sensing approach has been developed for monochromatic light. Broadband light leads to a radial 
smearing of the wavefront information as the speckle pattern scales with wavelength A. One way to overcome this 
problem is the use of an extremely chromatic optical system, which has a wavelength-dependent magnification 
such that the scaling with A due to diffraction is exactly compensated. Such optical systems have been described 
in the literature U^lliZI The broadband wavefront sensing will therefore require its own optical system, and we 
cannot directly use the science focal-plane for the wavefront sensing. 

The algorithm is currently limited to a single point source, which is often the case for exoplanet imaging, 
and can easily be implemented in laboratory applications. An application to extended objects is feasible, but it 
is not clear at this time that all wavefront modes can indeed be measured. 

Another requirement is that the aperture has an even symmetry. It therefore also works with typical 
amplitude-apodization pupil masks, but will not work with typical pupil phase-mask coronagraphs, which are 
even in their real part and odd in their imaginary part. It might be possible to modify our approach to phase-mask 
coronagraphs, which will be part of our future efforts. 

We have not optimized our approach in any ways regarding the different free parameters or the gains in 
the feedback loop. Such optimizations will be the topic of future investigations. And, finally, we have not yet 
considered the influence of amplitude aberrations and deviations between the expected and the true, perfect, 
PSF on our algorithm. The successful laboratory experiments performed by our teamP however, give us hope 
that these potential errors may be overcome. 

7. APPLICATIONS 

The original idea for this development was driven by the need for an efficient wavefront sensing and extreme 
adaptive optics control algorithm for high-contrast imaging that scales more gracefully than the classical number- 
of- actuators-squared (iV 2 ) law. Such an algorithm is indeed essential for future exoplanet imaging instruments 
for the next generation of extremely large telescopes. Apart from this originally intended application, our 
approach is useful for many other applications such as non-common path aberration measurements and slow 
shape corrections for future high-contrast imaging space missions. 

Indeed, our approach is ideally suited to measure small, non-common path aberrations in AO-assisted instru- 
ment. To determine them, one uses a monochromatic point source in a focus in front of the AO system. Such 
a light source is often already present to calibrate the adaptive optics itself. Our approach now provides for an 
iterative wavefront sensing in the focal plane using the deformable mirror to introduce the required diversity. 
Results from an actual test of this approach can be found in these proceedings!^ 
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